Reheating temperature from the CMB 
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In the recent paper by Mielczarek et al. (JCAP 1007 (2010) 004) an idea of the method which 
can be used to put some constraint for the reheating phase was proposed. Another method of 
constraining the reheating temperature has been recently studied by Martin and Ringeval (Phys. 
Rev. D 82 (2010) 023511). Both methods are based on observations of the cosmic microwave 
background (CMB) radiation. In this paper, we develop the idea introduced in this first article 
to put constraint on the reheating after the slow-roll inflation. We restrict our considerations to 
the case of a massive inflaton field. The method can be, however, easily extended to the different 
inflationary scenarios. As a main result, we derive an expression on the reheating temperature 
Trh. Surprisingly, the obtained equation is independent on the unknown number of relativistic 
degrees of freedom g* produced during the reheating. Based on this equation and the WMAP 7 
observations, we find Trh = 3.5 • 10 6 GeV, which is consistent with the current constraints. The 
relative uncertainty of the result is, however, very high and equal to ct(Trh)/Trh 
this uncertainty will be significantly reduced with future CMB experiments. 
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I. INTRODUCTION 

The reheating [J [2j is a hypothetical process in which 
the inflaton field [3] is converted into the standard model 
particles. The mechanism of reheating is usually assumed 
to be a parametric production of particles [4|. How- 
ever, the considerations are purely speculative due to the 
lack of any possible empirical verification of the reheating 
phase. Up to now, only some weak constrains on the ba- 
sic parameters of reheating are available. In particular, 
the reheating temperature Trh can be constrained from 
the both sides. From bottom the constraint is given by 
the big bang nucleosynthesis (BBN), namely Trh ^ 10 
MeV 0, H[. From the top, the constraint comes from 
the energy scale at the end of inflation Trh ;$ 10 16 GeV. 
Roughly 18 orders of magnitude remain to place the re- 
heating temperature somewhere between. Worse, there is 
no observational window available at these energy scales. 
Such a window exists however at the energies of infla- 
tion. It is because the perturbations created during the 
inflation can be studied by its impact on the cosmic mi- 
crowave background (CMB) radiation and subsequently 
by the large scale structures (LSS). The method of con- 
straining the reheating phase indirectly by the inflation- 
ary observational window was recently studied in Ref. 
• It was shown that it leads to the lower constraint on 
the reheating temperature Trh ^ 6 TeV. 

In this paper, we present an alternative method which 
can be used to fix at least some details of reheating based 
on observations of the cosmic microwave background 
(CMB) radiation. In comparison with the method pre- 
sented in Ref. 7], it will be possible not only to put a 
constraint on Trh, but just fix its value. The idea of the 
method was sketched in Ref. Q. It bases on the fact 
that the total increase of the scale factor from the ob- 
served part of inflation till now can be determined from 



the CMB. The number of e-foldings from the observed 
part of inflation till its end can be determined too. Based 
on this, the e-folding number from the end of inflation till 
the recombination can be found. As we show, this can be 
used to determinate the reheating temperature. For sim- 
plicity, we assume the slow-roll inflation (described by a 
massive inflaton field), which is in good agreement with 
the CMB observations. After inflation, the inflaton field 
undergoes coherent oscillations at the bottom of a po- 
tential well. The reheating takes place when the Hubble 
parameter H falls to the value of the inflaton decay rate 
r^. We assume that the reheating is instantaneous. Af- 
ter reheating, the standard radiation phase takes place. 
The evolution of radiation is assumed to be adiabatic. 
During the reheating, the effective number of relativis- 
tic species produced is given by g*. The decay rate of 
the inflaton field can be related with the remaining two 
parameters and Trh by the Friedmann equation as 
follows 
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where mp\ — 1.22 • 10 19 GeV. Therefore, only two from 
the parameters of reheating (r^,, Trh, <7*) are indepen- 
dent. In this paper, we show that the reheating temper- 
ature can be determined independently on the remaining 
two parameters. Up to now, the constraints on Trh were 
dependent on the value of g*. However, in the equation 
derived in this paper, the g* factors surprisingly cancel 
out. Having Trh, the decay rate can be expressed in 
terms of g* only. 

The considerations presented in this paper are re- 
stricted to the simplest setup in order to capture the 
essence of the method. However, extension to the differ- 
ent inflationary scenarios and to the more detailed mod- 
els of reheating can be done straightforwardly. 
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II. METHOD 



The main idea of the method can be understood by 
looking at Fig. [T] In this figure we schematically present 
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FIG. 1: Schematic evolution of the Hubble radius (blue line) 
for the standard cosmological scenario. The straight (red) line 
represents evolution of the physical length scale X(a), where 
X(a ) = Ao. 

evolution of the Hubble radius Rn := 1/H, together 
with the evolution of an arbitrary physical length scale 
A. The present value of this length scale is equal to Ao, 
what we call the pivot scale. The following values of the 
scale factor were distinguished: 

ao - the present value of the scale factor, we set ao = 1 
for convenience. 

ai - the scale factor at the end of the radiation era. 
Later, we set this to be a scale factor at the recombination 
which takes place soon after the end of the radiation era. 

a2 - the scale factor at which the instantaneous reheat- 
ing takes place (beginning of the radiation era) . 

as - the scale factor at the end of inflation. The infla- 
tion field starts to oscillate. 

a4 - the scale factor at which the length scale of 
the present value Ao crossed the Hubble radius during 
inflation. 

The total increase of the scale factor from 04 to ao will 
be of particular importance. We call it A to t, which can 
be expressed as follows 



Atot^J^Ai, where A; 



i=0 



(2) 



So, if we know durations of the four stages between the 
04 and ao, the A to t can be determined. This is, however, 
practically impossible to obtain because we do not know 
details of the intermediate periods as A 2 and Ai. Hope- 
fully, there is an alternative method to determinate A to t, 
which can be used to put constraints on A2 and Ai. This 
method bases on the observation of the CMB radiation. 
In particular, on the measurements of the scalar power 



spectrum. The form of this spectrum is parameterized 
by the function 



V s (k) = A s 



(3) 



The A s is an amplitude and n s is a spectral index of the 
scalar perturbations. The kg is some arbitrary fixed scale 
called pivot number. We can relate it to the pivot scale 
Ao introduced earlier by Ao = 2n/ko. In particular, the 
WMAP collaboration choice is k = 0.002 Mpc -1 (we 
also use this choice in this paper). For this value, the 
seven years of observations made by the WMAP satellite 
give the following values of the amplitude and spectral 
index of the scalar perturbations [9| 



A s = 2.441^-10- 



0.963 ±0.012. 



(4) 
(5) 



On the other hand, the well known prediction of the slow- 
roll inflation is 
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where e is the so-called slow-roll parameter equal to 
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For the considered massive slow-roll inflation n s = 1 — 4e. 

Let us now consider the power spectrum at the length 
scale Ao which corresponds to the pivot number fco- From 
observation, an amplitude of the scalar perturbations at 
this scale is equal to V a (ko) — A s . On the other hand, 
this amplitude is formed when k ~ aH. Therefore, for 
the mode ko we have S — A s . 

Cosmological evolution of the pivot scale Ao is given 

by 



A (a) = A — . 

ao 



(8) 



This relation is represented by the red line in Fig. [T] The 
value of A was equal to the Hubble radius at a\. Based 
on this, one can derive 
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where H is the value of the Hubble parameter when the 
A crossed the horizon during the inflation. In the second 
equality, we have used relation A (a) = ^r^j-, together 
with k ~ aH at the horizon crossing. Namely, A = 



and A(a4) = 7f^, where ao 
S = A s , so 



1. At the pivot scale, 
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Expressing the e from n s = 1 — 4e, we find 
1 mpi 



2 k 



(11) 



The essential conclusion derived from this equation is 
that: Based on the CMB observations, one can determi- 
nate the total increase of the scale factor from the ob- 
served moment of inflation till now. In principle, from 
the WMAP 7 observations we determinate 
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At„t = (8.0 ±1.5) • 10 



III. INFLATION AND REHEATING 



(12) 



The increase of a scale factor during the part of infla- 
tion from ct4 to <i3 is given by 



A, = 



,^0 



(13) 



where iV bs is the e-folding number, which can be ex- 
pressed as follows 
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We have used here V{4>) = ^(f> 2 and defined </> bs = 
0(04). In particular, based on the WMAP 7 data one 
can find iV b s = 54 ± 18. The uncertainty is high because 
of the strong sensitivity on the uncertainty of the spectral 
index n s . This will later propagate to the uncertainty of 
the reheating temperature. As we show in Sec. [V] the 
uncertainty of iV b s can be significantly reduced with the 
future CMB experiments. 

During the slow-roll inflation, evolution of the inflaton 
field cj> is well approximated by 



0(t) = 4>n 



771 mpi 



t. 



(15) 



From comparison with the numerical results, it can be 
seen that this approximation holds till the end of infla- 
tion, when <f> 0. Therefore, the kinetic term 
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is approximately constant during the inflation. This con- 
tribution to the total energy density is, however, domi- 
nated by the potential part during the slow-roll inflation. 
At the end of inflation, the contribution from the poten- 
tial part falls to zero (V(<j) = 0) = 0), and the kinetic 
term dominates. One can therefore estimate that, at the 
end of inflation, the energy density is given by 
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(17) 



A validity of this approximation was confirmed by the 
numerical computations. 

After inflation, the field starts to oscillate at the bot- 
tom of the potential well. During this evolution, the en- 
ergy density drops as in the matter dominated universe 
.10] 
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This evolution holds till when H ss and the re- 
heating takes place. Then, the energy density 



p(a 2 ) = g. 
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here g* — <?(Trh) is the number of ultrarelativistic de- 
grees of freedom generated during the reheating, where 
g{T) is defined as follows 
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In particular, for Glashow-Weinberg-Salam (GWS) 
model SU(2) L ® U Y (l) ® SU C (3), we have g = 106.75. 
Therefore one may expect that g* > 106.75 if the temper- 
ature of reheating is greater than the electroweak energy 
scale, T RH > 300 GeV. 
Based on (fTB")) we have 



A 2 = ^ = 



(12 
3 
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and applying (IT7|) and P^|) we derive 
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(22) 



This result will be useful in the subsequent section when 
deriving the expression on Trh- However, before we pro- 
ceed to this issue we can see what we already can say 
about the reheating temperature. Let us notice that the 
following condition p(a 2 ) < p{a^) must be fulfilled. En- 
ergy scale of reheating cannot be higher than energy at 
the end of inflation. In order to use this constraint one 
has to firstly determinate inflaton mass in Eq. p7|) . It 
can be done by noticing that the Friedmann equation 
reduces to 



H 
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in the slow-roll regime (e <C 1). Based on this and con- 
dition S = A Sl together with 77 s = 1 — 4e, one can derive 



= 777pii^37TA s (l - 77 s ). 



(24) 



Applying this expression to the WMAP 7 results we ob- 
tain 



777 = (1.4 ±0.5) • 10~ 6 777 P 1 

= (1.7 ±0.6) • 10 13 GeV. 



(25) 
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With use of this value, the condition piaz) < pips) re- 
duces to 



5* 1/4 T RH < 6.5 • 10 15 GeV 



(26) 



Based on this, part of the parameter space (3*,Trh) can 
be excluded. It was shown as the shadowed region above 
the thick line in Fig. [2] As we mentioned earlier, if 




300 



FIG. 2: Section of the reheating parameters space (Trh,5*). 
The shadowed region is excluded by Eq. [26] and g, > 106.75 
valid for Trh > 300 GeV. The thick line represents Trh = 



-1/4 



6.5 ■ 10 15 GeV. 



Trh > 300 GeV then > 106.75. This constraint ex- 
cludes another part of the parameter space. This was 
represented in Fig. [2] as the shadowed region constrained 
by the vertical line. Based on the above constraints, one 
can conclude that 



T RH < 2.0 • 10 15 GeV. 



(27) 



This corresponds to the inflationary bound on the reheat- 
ing temperature. 



IV. REHEATING TEMPERATURE 

After reheating, the Universe is filled by the relativistic 
plasma. Expansion of this relativistic gas is assumed to 
be adiabatic and the masses of particles are neglected. 
The adiabatic approximation is valid until the entropy 
transfer between the radiation and other components can 
be neglected. In turn, this second approximation is valid 
if the temperature is much higher than the masses of the 
particles. Then, dS = 0, which implies sa 3 = const, 
where the entropy density s of radiation is given by 



2ir 2 , 
45 



(28) 



Based on this one can derive expression on the increase 
of the scale factor from reheating till the recombination 
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We have T2 = Trh and T\ is equal to the recombination 
temperature T roc . During recombination gi = g 7 = 2 
and during reheating 52 = 9* , therefore 



Ax = 



rr.c 



(I) 
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Finally, increase of the scale factor from recombination 
till now is given by 



An = 1 + Z T , 



(31) 



where z rec is the recombination redshift which can be de- 
termined from the CMB observations. It is worth men- 
tioning that an intermediate stage other than recombi- 
nation can be used here. In particular, the equilibrium 
point (end of the radiation epoch, where p rac i = p ma t) 
can be chosen. The corresponding value of redshift can 
be also determined from the CMB observations. 

At this point, we have all required to find the expres- 
sion on Trh- We have found all Aj and A to t- Based on 
©, the following relation is fulfilled 



= A 3 A 2 AiA . 
Inserting (T3J), flSSJ, (J3DJ) and we obtain 
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(33) 



where we have used T roc = Tcmb(1 + z Tec ). The impor- 
tant observation is that g* factors cancel out. This is 
crucial, because the expression on Trh will be free from 
the dependence on the unknown g* parameter. With use 
of (TTTj) . (Ti"4"|) and the above equation can be rewrit- 
ten into the following form 



Trh 
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(34) 

This equation is a main result of this paper. Taking the 
constant parameters T CM b = 2.725 K = 2.348 • 10" 4 eV 
and fco = 0.002 Mpc -1 (and reexpressing units: Mpc -1 = 
6.39 • 10 -30 eV) one can rederive Eq. (|34| to the practical 
form 



Trh = 3.36 • 10 



-6M 



1 



A, 



■ exp 



6 



1 



GeV. (35) 



In Fig. [3] we show relation (|35|) as a function of the 
spectral index n s . We also mark the regions excluded 
from the inflationary constraint and the BBN constraint. 
For the data from the WMAP 7 observations, Eq. ([33)1 
leads to 



Trh = 3.5 • 10 6 GeV. 
The relative uncertainty of this result is 
ct(Trh) 



Trh 
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FIG. 3: The thick line represents Eq. [35] with A s = 2.441 ■ 
10 -9 . Dashed lines corresponds to n s — 0.963, taken from the 
WMAP 7 observations, what leads to T RH = 3.5 ■ 10 6 GeV. 
The shadowed regions are excluded by the inflationary and 
BBN constraints. 



Here a first order Taylor expansion was applied when 
calculating propagation of uncertainties: 



o-(Trh) 




a 2 (n s ) + 
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However, due to the strong (exponential) dependence of 
Xrh on n s , the applied linear approximation may turn 
out to be insufficient. Therefore, one can expect greater 
uncertainty of Trh than obtained here. Future studies 
need to address this issue. The high relative uncertainty 
(f3"T|) is mainly a result of the weakly determined value of 
Nobs, which is a function of n s . In the next section, we 
will examine how this uncertainty can be reduced with 
the future CMB experiments. 

As it was discussed in Refs. [n], [llj], if Trh ~ 
10 6 ~ 9 GeV, then it may be possible to measure Trh 
with the planned space-based laser interferometer ex- 
periments such as the Big Bang Observer. The value 
Trh = 3.5 • 10 6 GeV obtained here fulfills this condi- 
tion. Our prediction has therefore chance to be verified 
in future. 

Furthermore, based on Eq. (JTJ), one can express the 
inflaton decay rate in terms of <7* : 



1.7-10" 



GeV, 



(39) 



where the previously derived value of Trh was used. It 
is reasonable to expect that — a ■ m, which comes 
from the Heisenberg uncertainty relation. The a is a 
dimensionless parameter. With use of the inflaton mass 
found earlier wc find 



a 



10 



-19 



'9* 



(40) 



Based on this result, one can deduce that the inflaton 
decays into the very light particles comparing with its 
mass. This is also the reason why the reheating takes 
place at the relatively low energies. It becomes unsta- 
ble only when the sufficiently low energies are reached. 



However, the origin of this low vale of decay rate can- 
not be understood without a deeper understanding of the 
inflationary cosmology. 



V. FORECASTING 

As we have shown in the previous section, the value of 
Trh is strongly dependent on n s . Therefore, the method 
presented can be used effectively only if the value of n s 
is determined with high precision. The value of n s from 
the WMAP 7 observations is not determined sufficiently 
precise to obtain a strong prediction concerning the re- 
heating temperature. However, it may change if the new 
observational data will be available. In this section, we 
predict how the uncertainty on Trh will be reduced with 
the future CMB experiments. In particular, we consider 
the Planck satellite [l3| experiment which is currently on 
the stage of collecting data. We consider the ACTPol 
[HI ground-based experiment which is under construc- 
tion at present. We also consider the planned CMBPol 
[lj| satellite experiment. 

The uncertainty of Trh comes mainly from n s , there- 
fore, in the considerations we fix the value of A s . Follow- 
ing Ref. [jq . the expected uncertainties of n s from the 
mentioned CMB experiments are the following 



0.0031 Planck 
a{n s ) = { 0.0021 Planck+ ACTPol. 
0.0014 CMBPol 



(41) 



Based on this, let us first see the resulting uncertainties 
of the e-folding number A^bs- We find 



ct(AUs) 



Planck 
Planck+ACTPol. 
CMBPol 



(42) 



This significant reduction of the uncertainty of N ODS (with 
respect to the WMAP 7 results) will be crucial for deter- 
mining the reheating temperature. Based on (|35[) with 
(|¥T|) . we forecast 



Trh 




Planck 
Planck+ACTPol. 
CMBPol 



(43) 



Here, the values of the parameters n s and A s were set to 
be those obtained from the WMAP 7 observations. 

In order to have a (Trh) /Trh smaller than unity, the 
uncertainty of n s should be reduced by 2 orders of mag- 
nitude with respect to the WMAP 7 results. At present, 
there is however no experiment planned to reach such 
sensitivity. The uncertainty may be nevertheless addi- 
tionally reduced by combining data from the different 
available experiments. This is possible because of the 
angular scale dependent sensitivity of the CMB experi- 
ments. In particular, the ground-based experiments can 
provide much better data of the CMB polarization at the 
small angular scales (high multipoles) than the space- 
based experiments can. 
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VI. SUMMARY 

In this paper, we have developed a new method of con- 
straining the reheating phase after inflation. The method 
bases on the observations of the cosmic microwave back- 
ground radiation. In particular, the fact that the total 
increase of the scale factor from the observed part of infla- 
tion till now can be determined is used. Based on this, we 
have found the expression on the reheating temperature. 
The expression is free from the dependence on the un- 
known parameter. With use of the WMAP 7 results, 
we have determined Trh = 3.5 • 10 3 TeV. The relative 
uncertainty of this result is equal to ct(Trh)/7rh ~ 53. 
This high uncertainty can be, however, significantly re- 
duced with the future CMB data. One can expect 
the reheating temperature to be quite precisely deter- 
mined (reaching ct(Tr,h)/Trh = 0(1)) within the present 
decade. 

The value of reheating temperature determined in this 



paper is consistent with the known bounds, in particular, 
with the lower bound Trh ^ G TeV recently found in Ref. 
7]. It is also interesting to note that within the super- 
symmetric extension of the standard model, the upper 
bound on the reheating temperature exists Trh ^ 10 4 
TeV (see Refs. [1343 )■ Our result is also in agreement 
within this condition. Finally, it is worth mentioning 
that the low value of the reheating temperature, as de- 
termined here, can have interesting implications on the 
phenomenology of primordial black holes [20j . 
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